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ABSTRACT 

A  one-dimensional  mixed-layer  model  of  ocean  thermal 
structure  (Camp,  19  76)  was  first  modified  to  determine  the 
effects  of  a  salinity  profile  on  density  structure  and 
secondly  to  observe  the  added  effects  of  surface  salinity 
flux.   The  model  was  tested  in  August  and  November-December 
1974  with  (1)  hydrocast  data  at  Ocean  Stations  HOTEL  and 
PAPA  and  (2)  salinity  profiles  statistically  derived  from 
historical  data  for  six  ocean  stations  and  adjusted  to 
correspond  to  an  observed  BT.  Inclusion  of  salinity  struc- 
ture did  not  significantly  affect  mixed-layer  depth  or 
temperature  predictions  during  the  summer  period.   During 
the  winter  period,  salinity  structure  tended  to  inhibit 
deepening  yielding  a  slightly  warmer,  shallower  mixed 
layer.   Surface  salinity  flux  altered  significantly  the 
thickness  of  the  isothermal  layer,  with  decreases  in  thick- 
ness during  precipitation  (E-P<0)  and  a  tendency  for  in- 
creased thickness  during  some  periods  of  net  downward  heat 
flux  when  E-P>0 . 
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I.   INTRODUCTION 

It  has  been  established  that  the  upper  layer  of  the 
ocean  responds  to  atmospheric  forcing  on  the  synoptic  and 
diurnal  time  scales  [Elsberry  and  Camp  (19  78),  Garwood 
(19  77)].   The  success  of  those  models  which  incorporate 
only  mechanical  mixing,  convective  overturning,  and  vertical 
diffusion  in  predicting  the  time  evolution  of  the  mixed 
layer  (Camp  and  Elsberry,  19  78)  has  shown  the  short  term 
ocean  boundary  layer  response  to  be  predominately  one-dimen- 
sional in  many  ocean  regimes .   Although  the  general  model 
[see  Camp  (19  76)]  contains  both  temperature  and  salinity 
effects,  Camp  and  Elsberry  (19  78)  considered  density  as  a 
function  of  temperature  only,  and  assumed  buoyancy  flux  to 
the  atmosphere  to  be  attributable  to  the  heat  flux  alone. 
Johnson  (19  77)  followed  with  an  evaluation  of  the  model's 
ability  to  forecast  ocean  thermal  structure  based  on  atmos- 
pheric forcing  derived  from  synoptic-scale  data.   Again, 
due  to  lack  of  salinity  data,  Johnson  neglected  salinity 
effects  and  therefore  had  to  adjust  initial  temperature  pro- 
files to  remove  positive  gradients.   It  is  the  purpose  of 
this  research  to  study  salinity  effects  on  the  upper  ocean 
thermal  structure  using  the  Camp  (19  76)  model  by  making  com- 
parisons with  the  results  of  Johnson  (1977). 

Miller  (1976)  included  salinity  in  the  mixed  layer  model 
of  Kraus  and  Turner  (1967).   Using  BOMEX  Period  III  upper 
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ocean  soundings,  he  showed,  for  periods  up  to  10  days,  the 
effects  of  salinity  on  mixed-layer  temperatures  and  depths 
and  the  effects  of  precipitation  on  surface  salinity  and 
temperature  in  the  tropical  Atlantic.   He  concluded  that 
the  most  significant  effect  of  inclusion  of  salinity  was 
the  reduced  cooling  rate  of  the  mixed  layer  due  to  reduced 
entrainment  of  water  from  the  pycnocline  into  the  base  of 
the  layer.   This  is  a  result  of  a  larger  density  jump  at 
the  base  of  the  layer  when  a  salinity  jump  corresponds  to 
a  temperature  jump  than  would  be  the  case  if  density  were 
a  function  of  temperature  alone.   This  type  of  density 
structure  can  also  lead  to  a  positive  temperature  gradient 
(temperature  increasing  with  depth)  immediately  below  the 
mixed  layer  (i.e.  the  mixed-layer  temperature  is  cooler 
than  that  of  the  underlying  water).   Miller  also  concluded 
that  the  effects  of  precipitation  were  greatest  under  condi- 
tions of  light  winds  and  heavy  rainfall  during  which  a 
shallow  stable  layer  is  produced  at  the  ocean  surface. 
Considering  Miller's  findings,  it  is  apparent  that 
salinity  is  important  in  determining  the  density  structure 
when  applying  a  mixed-layer  model  to  the  entire  ocean.   Fur- 
thermore, although  the  relative  importance  of  salinity 
effects  on  the  short  term  evolution  of  the  density  profile 
may  not  be  significant,  one  would  anticipate  that  a  long 
term  imbalance  between  evaporation  and  precipitation  contri- 
butes to  the  buoyancy  flux  and  could  lead  to  significant 
density  changes.   Additionally,  sound  speed  in  the  ocean 
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varies  by  1.2  to  1.3  ms~   (  /oo)    and  salinity  variations 
in  the  near-surface  layers  can  thus  have  significant  ef- 
fects on  the  sound-speed  profile.   These  effects  would  be 
greatest  in  regions  of  varying  surface  salinity  flux  and 
advection  (e.g.,  doldrums,  subpolar  regions). 

The  importance  of  salinity  in  the  model  and  the  paucity 
of  concurrent  temperature  and  salinity  data,  even  at  ocean 
weather  stations,  raise  the  question  of  how  to  incorporate 
salinity  effects  into  the  initial  conditions  and  to  then 
verify  the  results  when  the  only  available  information  is 
a  temperature  profile  from  a  bathythermograph.   It  was  the 
objective  of  this  study  to  test  the  salinity  contribution 
to  the  density  structure  observed  based  on  two  different 
approaches:   (1)  utilization  of  available  hydrocasts  when 
both  temperature  and  salinity  profiles  are  available,  and 
(2)  utilization  of  bathythermographs  and  typical  tempera- 
ture/salinity profiles  derived  from  the  historical  data 
files  and  adjusted  to  be  compatible  with  an  observed  BT. 
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II.   MODEL  MODIFICATIONS 

Camp  (19  76)  and  Miller  (19  76)  presented  the  governing 
equations  for  the  vertical  redistribution  of  heat  and  turbu- 
lent kinetic  energy  in  the  upper  ocean.   These  equations,  as 
applied  to  a  well-mixed  layer,  were  simplified  by  appropri- 
ate assumptions  and  closure  was  achieved  by  parameterization 
of  the  dissipation  and  shear  production  terms  of  the  turbu- 
lent kinetic  energy  budget.   At  the  surface,  the  boundary 
conditions  were  specified  in  terms  of  the  wind  stress ,  inso- 
lation, sensible  and  latent  heat  fluxes,  back  radiation,  and 
evaporation  minus  precipitation.   The  boundary  conditions  at 
the  base  of  the  layer  are  then  computed  by  the  model  in 
terms  of  the  fluxes  of  momentum,  heat,  and  salt.   In  develop- 
ing the  numerical  scheme  from  the  expressions,  Camp  neglected 
salinity  effects  because  of  a  lack  of  salinity  data.   For 
this  study,  the  model  was  modified  to  account  for  the  contri- 
butions of  salinity  to  the  density  structure  and  the 
buoyancy  flux. 

A  linear  equation  of  state  was  assumed,  which,  when  ex- 
pressed in  terms  of  buoyancy,  B,  is 

B  =  agT  -  BgS  (2-1) 

where 

a  =  -  —   3T   and   6  =  —  TS 
o  Ho 
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_4  o   -1 
are  taken  as  constant;  a  =  2 . 5  x  10   (  C)    and  g  = 

7.5  x  10"I|(O/oo)~:L  after  Camp  (1976).   The  assumption  that 
3  is  constant  is  quite  good.   However,  density  is  not  a 
linear  function  of  temperature  and  a  can  be  taken  as  con- 
stant only  over  a  limited  range  of  temperatures.   The  den- 
sity profile  is  represented  in  the  model  as  temperature  and 
salinity  pairs  to  a  depth  of  -NAZ  by  N  isothermal  and  iso- 
haline  slabs  of  AZ(=  2.5  m). 

The  potential  energy  (PE)  per  unit  is  then  expressed  as 

o 


/ 


PE  =  -pQg  J    (aT  -  gS)  ZdZ  (2-2) 

-D 

where  D  is  a  fixed  depth  greater  than  the  maximum  penetra- 
tion of  the  vertical  turbulent  processes .   Assuming  advec- 
tive  processes  and  the  variation  due  to  compressibility  are 
negligible,  the  change  in  potential  energy  resulting  from 
mixing  the  layer  from  a  depth  -NAZ  to  a  depth  -(N+DAZ  is 


APE(N)  =  h    pogN(AZ)2[a(TN-TN+1)-B(SN-SN+1)]  .         (2-3) 


If  a(T„-T^  .  )-g (S„-S„  ,  )  >  0,  the  water  column  is  stable 
and  APE  (N)  is  the  amount  of  turbulent  kinetic  energy  re- 
quired  to  mix  the  layer  to  a  depth  -(N  +  DAZ.   If  a(T„-T    ) 
-  g(SM-S„+,)  <  0,  the  column  is  unstable  and  will  overturn, 
generating  turbulent  kinetic  energy  by  free  convection  (re- 
presented as  APE  (N)). 

c 
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In  each  one  hour  time  step  the  model  first  accounts 
for  surface  fluxes ,  and  the  resulting  temperature  profile 
and  surface  (N=l)  salinity  are  determined  by 


T  (t+At)  =  T1(t)+-^|-[-QT(0,t*)+Q  (0,t*)-Q  (-AZ,t*)] 

o  p 

TM(t+At)  =  TM(t)+  -^-[Q  (-NAZ,t*)-Q  (-  (N  +  l )  AZ  ,  t  *  )  ]       (2-4) 
N  N     p  c   xs  s 

o  p 

S1(t+At)  =  S1(t)+  S1(t)  ||[E(t*)  -  P(t*)] 


where  Q   is  solar  radiation,  QT  is  the  sum  of  latent,  sen- 
sible,  and  back  radiation,  E  is  the  evaporation  rate,  and 
P  is  the  precipitation  rate  during  the  interval  t*  =  At. 
Insolation  is  distributed  vertically  in  the  water  column 
by  absorbing  50%  in  the  first  meter  and  absorbing  the  remain- 
ing 50%  as  exp(-yZ).   The  extinction  coefficient,  Y,  was 
assumed  constant  at  0.002  cm    after  Johnson  (1977). 

If  the  resulting  buoyancy  profile  is  unstable,  the  water 
column  is  mixed  according  to 


TN>11)  =  (NT^l}  +  T^)/(N  +  1) 


S(i+1)  =  (NS(i)  +  S(i+1)/(N+1) 
N  +  l         1       N  +  l    UNJ  +  1; 


(2-5) 


until  the  profile  is  stable.   The  temperature  and  salinity 
profiles  are  then  isothermal  and  isohaline  to  a  depth 
-(N  +  DAZ.   The  potential  energy  generated  by  this  free  con- 
vection is  found  from 


°  NN-1 

APE  =-p  gl   (ct6T-B<5S)ZdZ«  /^ 
r,  N=l 


APE  (N)   (2-6) 
c 


-D 
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where  6T  and  5S  are  the  changes  in  temperature  and  salinity 
resulting  from  mixing  and  the  depth  of  free  convection  is 
(NN)CAZ)  <  D. 

If  the  resulting  buoyancy  profile  is  stable,  additional 
kinetic  energy  must  be  expended  to  further  mix  the  water 
column.   The  total  amount  of  energy,  E„,  available  for  mix- 
ing the  first  N  levels  with  the  N+l  level  is 

E,=E+E-E  (2-7) 

T    m     c     p 

where  E   is  the  mechanically  generated  turbulent  kinetic 

m  J    to 

energy,  E  is  the  convectively  generated  turbulent  kinetic 
energy,  and  E  is  the  amount  of  energy  previously  expended 
to  mix  the  water  column  to  depth  -NAZ. 

E   is  calculated  by 

m  J 

E  (N)  =  [p  w...3  exp(-MAZ/H)  ]  At  (2-8) 

m        o  -     r 

with  w.,.  =  average  friction  velocity  in  the  ocean  over  the 
time  interval  At  and  H  =  scale  depth  of  50  m.   This  depth- 
dependent  formulation  follows  the  approach  of  Elsberry, 
Fraim,  and  Trapnell  (1976).   Setting  a  scale  depth  of  H=°° 
would  correspond  to  the  original  Kraus  and  Turner  (1967) 
model  in  which  the  net  generation  minus  dissipation  was  in- 
dependent of  the  depth  of  the  layer. 

E   is  considered  to  be  the  fraction  of  APE   available 
c  c 

for  entrainment  and  is  calculated  as 

NN-1 

E  (N)  =  -R  y   APE  (M)  (2-9) 

c  /  j  c 

N=l 
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where  R  =  0.15  following  Camp  (1976)  and  Johnson  (1977), 

and  represents  the  percentage  of  convectively  generated 

energy  not  dissipated  by  viscous  forces.   E   is  calculated 

from 

N-l 

E  (N)  =  Y^    APE  (i)  (2-10) 

p       L-d  m 

i  =  NN 

If  Em  >  APE   for  N  levels,  then  there  is  energy  avail- 
T      m  ° 

able  to  mix  to  level  N+l.   If  E^  <  APE  ,  the  level  is  par- 

T       m 

tially  mixed  following  Thompson  (1976)  and  Camp  (1976). 
Thus ,  let 

NET(N) 

a  =  (N+DAPE  (N) 
m 

Q  =  [aQN+1  -  (N-a)  QN]/N  (2-11) 

Q'  =  aQN  +  (1-a)  QN+1 

and  set  Q.  =  Q  for  i=l,2,...,H  and  Qvr  +  -,  =  Q'  where  Q  is  tem- 
perature and  salinity.   Note  that  the  mixing  algorithms  re- 
quire that  the  isothermal  and  isohaline  layers  be  the  same 
depth.   The  mixed-layer  depth  is  determined  from  the  tempera- 
ture profile  only  by  defining  it  as  the  depth  at  which  the 
temperature  is  0.2  C  less  than  T, ,  where  T   has  been  assumed 
to  be  equal  to  T  ,  the  sea-surface  temperature.   Light  wind 
conditions,  strong  heating,  or  heavy  precipitation  could  re- 
sult in  some  difference  between  T,  and  T  .   The  pycnocline 

1       o        ^J 

is  assumed  to  be  weakly  diffusive  and  the  change  of  tempera- 
ture and  salinity  below  the  well-mixed  layer  is  specified  by 
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II  -  a   32t 

,  (2-12) 

3S  .  A  jTS 

dt    '      v  8Z2 


where  the  vertical  eddy  diffusion  coefficient,  A  ,  is  taken 

2     -1 
to  be  0.5  cm  -sec    after  Haney  and  Davis  (1976).   This 

value  is  consistent  with  many  observations.   It  is  much 
greater  than  the  values  of  molecular  conductivity  and  dif- 
fusion.  The  processes  involved  are  not  well  understood  but 
may  be  due  to  breaking  of  internal  waves,  double  diffusion, 
and/or  lateral  advection  on  a  small  scale. 
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III.   DATA 

This  study  used  atmospheric  and  bathythermograph  data 
extracted  from  the  Fleet  Numerical  Weather  Central  (FNWC) 
historical  data  file  by  Johnson  (1977).   Additionally,  the 
12-hour  accumulated  precipitation  field  produced  by  FNWC's 
primitive  equation  model  was  extracted.   Oceanographic  sta- 
tion data  for  Ocean  Station  HOTEL  were  obtained  from  the 
National  Ocean  Data  Center  (NODC)  and  for  Ocean  Station  PAPA 
from  the  Institute  of  Ocean  Sciences  (1975).   The  data  at 
OS  HOTEL  consisted  of  18  Nansen  cast  observations  and  the 
data  at  OS  PAPA  consisted  of  three  hydrographic  observations 
and  11  STP  (salinity-temperature-pressure)  observations  for 
the  time  frames  indicated  in  Table  I.   Additionally,  daily 
surface  salinity  samples  were  available  at  OS  PAPA  for  the 
period  06  December  to  31  December  1974. 

To  take  advantage  of  the  more  numerous  bathythermograph 
(BT)  observations,  it  was  proposed  to  use  some  type  of 
historical  salinity  profile  in  conjunction  with  a  BT  to 
initialize  the  model.   One  possibility  was  to  use  a  profile 
generated  from  climatological  means  at  different  depths. 
This  proved  unfeasible  in  that  the  profiles  contained  very 
little  detailed  information  on  the  salinity  structure  in 
the  upper  levels  of  the  ocean.   Another  possibility,  using 
standard  temperature-salinity  relationships  for  a  water  mass, 
was  rejected  since  these  relationships  fail  in  the  boundary 


20 


layer  because  of  changes  caused  by  the  surface  fluxes.   The 
third  possibility  was  to  utilize  the  Hydroclimatological 
Data  Base  Retrieval  (HYDAT)  system  installed  at  FNWC  (Ocean 
Data  Systems,  Inc.,  1977).   HYDAT  retrieves  hydroclimatologi- 
cal data  from  the  historical  and  synoptic  data  bases.   FNWC ' s 
historical  files  contain  expendable  bathythermograph  and 
Nansen  cast  data.   The  synoptic  file  contains  the  most  re- 
cent 72  hours  of  BT  reports.   HYDAT  then  statistically  ana- 
lyses these  data  as  described  in  Appendix  A  and  selects  a 
temperature  profile  that  is  considered  typical  for  the  time 
period  (>  1  month)  specified. 

If  the  typical  profile  is  a  Nansen  cast,  the  salinity 
profile  is  output  with  the  temperature  profile.   If  the 
typical  profile  is  a  bathythermograph,  then  no  salinity  pro- 
file is  given.   In  this  case,  it  may  still  be  possible  to 
obtain  a  salinity  profile  by  searching  manually  through  the 
HYDAT  output  for  a  Nansen  cast  sounding  whose  sea-surface 
temperature,  deep  temperature,  and  heat  content  come  closest 
to  fitting  the  median  values.   For  this  study,  the  parameters 
were  weighted  such  that  the  fit  of  sea-surface  temperature 
and  heat  content  were  considered  more  important  than  the 
fit  of  deep  temperature.   In  some  locations,  Nansen  casts 
were  sparse  and  their  parameters  deviated  from  the  median 
values.   Thus,  a  minimum  requirement  for  considering  a  Man- 
sen  cast  as  a  typical  profile  was  that  the  three  parameters 
must  fall  within  the  range  between  the  first  and  third 
quartiles .   Proceeding  in  this  manner,  salinity  profiles 
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were  obtained  for  all  six  ocean  weather  stations  in  December 
and  for  three  ocean  weather  stations  in  August  (Table  I). 
These  salinity  profiles  were  assumed  to  be  typical  for  these 
ocean  station  locations  and  time  of  year.   More  ideally,  one 
would  like  to  have  a  data  bank  of  hydrocast  data  and  an  ex- 
tracting system  which  would  take  a  bathythermograph  and 
search  the  historical  data  for  the  closest  fitting  tempera- 
ture profile  and  then  use  the  corresponding  historical 
salinity  profile  as  the  typical  profile. 

Once  the  typical  profiles  were  obtained,  the  question 
arose  as  to  how  to  utilize  them.   The  initial  observed  tem- 
perature profile  likely  has  structural  details  that  should 
be  reflected  in  the  initial  salinity  profile.   These  details 
include  the  depth  of  the  well-mixed  layer  and  positive  tem- 
perature gradients.   Thus  the  typical  salinity  profile  should 
be  mixed  as  deep  as  or  deeper  than  the  BT  and  the  salinity 
gradient  should  be.  such  that  positive  temperature  gradients 
do  not  lead  to  an  unstable  density  profile.   The  typical  pro- 
files for  a  month  were  therefore  treated  in  two  ways  to 
generate  a  salinity  profile  that  would  correspond  to  any 
specified  BT  for  that  month. 

The  first  method  was  to  combine  the  typical  temperature 
and  salinity  profiles,  T„(Z)  and  S™(Z)  into  a  typical  buoy- 
ancy profile  BT(Z),  by 

BT(Z)  =  agTT(Z)  -  BgST(Z)  (3-1) 
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and  the  buoyancy  profile  was  checked  for  stability.   Then 
the  mixed-layer  depth,  the  depth  where  the  temperature  is 
0.2°C  less  than  the  surface  temperature,  of  the  typical  tem- 
perature profile  and  of  the  bathythermograph  were  determined. 
If  the  layer  depth  of  the  T^CZ)  profile  was  less  than  that 
of  the  bathythermograph,  the  buoyancy  profile  was  mixed 
according  to  equations  (2-5)  to  the  layer  depth  of  the  BT. 
This  mixing  is  necessary  because  of  the  assumption  in  the 
model  that  the  isothermal  and  isohaline  layers  are  mixed  to 
the  same  depth.   If  the  mixing  were  not  done,  the  model 
would  use  a  certain  amount  of  input  energy  in  mixing  the 
salinity  profile  to  the  level  of  the  temperature  profile, 
whereas,  in  reality,  it  is  likely  that  the  actual  salinity 
profile  was  already  mixed  to  the  deeper  level.   If  the  layer 
depth  of  the  T„(Z)  profile  was  greater  than  that  of  the  BT, 
no  mixing  was  done.   If  a  positive  temperature  gradient  was 
encountered  in  the  BT  before  the  layer  depth  was  reached, 
the  BT  was  considered  to  have  a  mixed-layer  depth  equal  to 
zero.   Finally,  with  any  given  bathythermograph  observation, 
T  (Z),  and  the  assumed  buoyancy  profile  Bm(Z),  a  correspond- 
ing salinity  profile,  S  ,  was  found  from 

agT  (Z)  -  BT(Z) 

S  (Z)  =  - .  (3-2) 

3g 

Therefore  this  procedure  generates  a  salinity  profile  from 
an  observed  profile  consistent  with  the  typical  buoyancy 
profile,  (3-1).   Since  this  buoyancy  profile  is  based  on 
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the  typical  temperature  and  salinity  profiles,  it  is  repre- 
sentative of  the  conditions  at  that  location  during  the 
specified  month. 

The  second  method  was  to  take  the  typical  salinity  pro- 
file and  match  it  directly  to  the  bathythermograph  observa- 
tion.  The  advantage  of  the  typical  salinity  profile  over  a 
climatological  profile  is  that  the  typical  profile  is  an 
actual  observation  and  therefore  has  detailed  structural 
features  that  the  climatological  profile  is  lacking.   The 
mixed-layer  depth  of  the  BT  was  determined  as  above  and  the 
typical  salinity  profile  was  required  to  be  mixed  to  the  same 
or  greater  depth.   The  resulting  profiles  gave  rise  to  a 
buoyancy  profile  which  was  then  required  to  be  stable: 
B(Z-AZ)  >_  B(Z).   If  the  profile  was  unstable  at  any  level, 
Z,  the  salinity  profile  was  adjusted  by 

<x[T  (Z)-T  (Z-AZ)] 
ST(Z)  =  ST(Z-AZ)  +  - g-2 (3-3) 

which  assures  neutral  stability.   A  new  B(Z)  was  calculated 
and  the  procedure  continued  until  the  entire  profile  was 
stable . 

The  salinity  profile  resulting  from  either  of  these 
methods  permits  positive  temperature  gradients  that  would 
appear  unstable  when  considered  alone.   Including  salinity 
therefore  allows  the  model  to  be  used  over  a  much  wider  range 
of  oceanic  conditions  without  biasing  potential  energy  as 
in  Johnson  (1977). 
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IV.   DISCUSSION  OF  RESULTS 

Initially,  the  salinity  flux  at  the  surface  of  the  ocean 
was  neglected.   This  is  potentially  accurate  only  if  there 
is  a  balance  of  evaporation  and  precipitation  (E-P)  over 
the  time  period  of  model  integration.   The  validity  of  this 
assumption  as  well  as  the  assumption  of  one-dimensionality 
can  be  checked  by  examining  Figs.  1-3.   In  Fig.  1,  the 
variability  of  surface  salinity  at  OS  PAPA  between  the 
first  STP  cast  (1724  GMT  2  3  November  19  74)  and  the  last 
(1730  GMT  06  December  1974)  is  quite  small  and  one  can 
safely  neglect  the  advection  and  surface  flux  of  salinity. 
However,  the  variations  of  temperature  and  salinity  in  the 
pycnocline  appear  to  be  characteristic  of  advective  affects 
associated  with  tidal  period  and  internal  wave  motion. 
Also  notice  that  the  positive  temperature  gradients  in  the 
layer  80-150  m  are  excellent  examples  of  the  need  for  cor- 
responding salinity  profiles  when  vertical  mixing  is  being 
modeled.   Variations  at  OS  HOTEL  between  1430  GMT  2  8 
November  and  1442  GMT  29  December  1974  (Fig.  2)  are  on  the 
order  of  one-half  part  per  thousand  and  generally  confined 
to  the  upper  50  m.   Notice  that  there  is  more  variability 
in  the  pycnocline  here  than  at  OS  PAPA.   At  OS  HOTEL  in 
August  (Fig.  3),  salinity  variations  are  quite  extreme  even 
within  individual  soundings.   The  variations  cannot  be 
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accounted  for  by  realistic  surface  salinity  fluxes  and 
therefore  are  likely  due  to  advective  effects  and/or  small 
changes  in  ship  position  and  the  proximity  to  the  Gulf 
Stream. 

A.   HYDROCAST  VS.  BT  INITIALIZATION 

The  first  model  run  shown  using  hydrocast  data  is  at  OS 
PAPA  during  the  period  from  1800  GMT  2  3  November  to  1800 
GMT  06  December  1974.   Wind  stress  and  total  heat  flux, 
Figs.  4  and  5,  respectively,  for  this  period  were  extracted 
from  the  FNWC  fields.   The  total  heat  flux  is  the  sum  of 
solar  radiation,  sensible  and  latent  heat  fluxes,  and  back 
radiation  and  is  generally  positive  for  this  period,  indi- 
cating an  upward  heat  flux  at  the  sea  surface.   Under  these 
conditions ,  one  would  expect  a  deepening  and  cooling  of  the 
mixed  layer.   The  initialization  and  verification  profiles 
from  STP  data  are  shown  in  Fig.  6  and  the  predicted  time 
evolution  of  surface  salinity,  surface  temperature  and  mixed- 
layer  depth  are  shown  in  Figs.  7,  8, and  9.   The  correspond- 
ing initial  temperature  profile  (with  positive  gradients 
removed  to  ensure  stability)  derived  from  a  BT  and  the  pre- 
dicted changes  in  mixed-layer  temperature  and  depth  are 
depicted  in  Figs.  10,  11,  and  12. 

Comparison  of  Figs.  8  and  11  shows  a  similar  cooling 
trend  and  good  agreement  with  the  data.   Notice  that  the 
initial  SST  from  the  STP  cast  was  slightly  warmer  than  the 
SST  from  the  BT.   Similarly,  in  Figs.  9  and  12,  one  can  see 
that  the  STP  initialization  predicts  a  consistently 
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shallower  MLD  as  compared  to  the  BT  initialization  but  each 
forecast  agrees  well  if  verified  against  similar  type  of 
data  (see  Figs.  6  and  10).   This  discrepancy  in  observed 
thermal  structure  is  seen  to  be  fairly  consistent  in  that 
the  sea-surface  temperature  averaged  0.2  6   C  warmer  and  the 
mixed-layer  depth  averaged  11.4  m  shallower  in  the  STP  cast 
data  than  the  BT  data  (Table  II).   Taking  this  into  account, 
the  model  performance  is  very  similar  throughout  the  period 
as  the  STP  initialization  resulted  in  a  cooling  of  the 
mixed  layer  by  0 . 8   C  and  a  deepening  of  9.8  m  whereas  the 
BT  initialization  resulted  in  a  cooling  of  0.78   C  and  a 
deepening  of  7.1  m.   Note  that  the  slab  depth  (depth  of  the 
isothermal  layer)  in  Figs.  9  and  12  retreats  and  deepens  in 
a  similar  manner  in  response  to  weak  stress  and  downward 
heat  flux  as  described  by  Johnson  (19  77).   This  would  be  an 
expected  feature  since  salinity,  in  the  absence  of  sources 
and  sinks  ,  would  have  no  effect  on  the  shoaling  of  the  iso- 
thermal  layer. 

The  major  contribution  of  salinity  in  this  case  is  the 
stabilizing  effect  on  the  positive  temperature  gradient 
below  the  thermocline.   The  salinity  gradient  that  allows 
this  to  occur  increases  the  density  jump  at  the  base  of  the 
mixed-layer  depth  and  tends  to  inhibit  deepening  of  the 
layer.   This  case  further  points  out  the  importance  of  the 
initial  profile  in  the  subsequent  model  predictions,  as  was 
also  shown  by  Johnson  (1977). 

The  next  case  considers  a  31-day  period  (2  8  November  to 
29  December)  at  OS  HOTEL.   The  synoptic  variability  of  stress 
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(Fig.  13)  and  the  corresponding  upward  heat  flux  (Fig.  14) 
again  indicate  a  deepening  regime.   The  initial  temperature 
profile  from  the  Nansen  cast  in  Fig.  15  and  from  the  BT  in 
Fig.  19  are  essentially  identical.   The  initial  salinity  pro- 
file (Fig.  15)  shows  a  relative  maximum  at  approximately 
100  m,  which  should  act  as  a  barrier  to  deepening.   However, 
below  this  level,  salinity  decreases  with  depth  and  would 
contribute  to  vertical  overturning  if  not  stabilized  by  the 
temperature  profile.   A  gradual  increase  of  surface  salinity 
is  predicted  (Fig.  16)  as  more  saline  water  is  entrained  at 
the  base  of  the  layer.   This  variation  appears  to  be  smaller 
than  the  natural  variability  in  the  salinity  measurements. 

Almost  identical  cooling  rates  of  the  surface  of  3.7   C 
during  the  period  were  forecast  by  hydrocast  initialization 
(Fig.  17)  and  BT  initialization  (Fig.  20).   Qualitatively 
this  follows  the  trend  indicated  by  the  observations ,  but 
forecasts  sea-surface  temperatures  that  are  on  the  order  of 
1.5  to  2.0   C  too  cold.   Likewise,  a  very  similar  trend  in 
mixed-layer  deepening  is  forecast  by  the  hydrocast  initiali- 
zation (Fig.  18)  and  by  the  BT  initialization  (Fig.  21). 
Both  forecasts  tend  to  indicate  too  deep  a  MLD  in  the  first 
14  days  as  compared  to  observations  but  then  tend  to  corre- 
spond to  the  data  in  the  later  half  of  the  period.   As  men- 
tioned previously,  the  slab  depth  responds  identically  in 
both  cases.   At  the  end  of  the  period,  the  Nansen  cast  ini- 
tialization is  approximately  7  m  shallower  than  the  BT 
initialization.   This  difference  is  due  solely  to  the 
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salinity  structure  because  the  same  forcing  is  used  in  each 
case  and  the  initial  temperature  profiles  were  nearly 
identical. 

It  may  seem  inconsistent  that  both  mixed  layers  would 
have  the  same  temperatures,  since  the  deeper  layer  would 
have  entrained  more  cold  water  from  below.   However,  given 
the  relatively  weak  temperature  gradient  and  a  difference 
in  depth  of  only  7  m,  the  extra  water  entrained  in  the 
deeper  layer  has  an  insignificant  cooling  effect  on  a  water 
column  that  is  on  the  order  of  100  m  deep. 

B.   TYPICAL  PROFILE  INITIALIZATION 

The  lack  of  adequate  temperature-salinity  observations 
over  large  parts  of  the  ocean  was  motivation  for  attempting 
to  use  the  salinity  profiles  derived  in  Section  III.   Both 
procedures  were  tested  in  order  to  determine  which,  if  either, 
was  more  useful. 

The  typical  buoyancy  profile  method  for  deriving  a 
typical  salinity  profile,  Equations  (3-1)  and  (3-2)  was 
tested  first.   This  method,  while  providing  initial  profiles 
that  were  stable,  gave  salinity  profiles  that  were  not  in 
general  representative  of  those  observed  in  the  ocean.   That 
is,  unrealistic  salinity  variations  were  created  below  the 
mixed  layer  due  to  the  constraint  that  buoyancy  be  conserved. 
This  is  shown  schematically  in  Fig.  22.   The  typical  buoyancy 
profile  in  Fig.  22a  is  mixed  from  Z,  to  Z~  to  correspond  with 
the  MLD  of  the  BT.   The  salinity  profile  derived  from  the 
typical  buoyancy  profile  is  shown  in  Fig.  22b.   In  the  mixed 
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layer,  buoyancy  is  constant,  thus  the  derived  salinity  will 
vary  according  to  the  temperature  by  Equation  (3-2).   At 
the  base  of  the  mixed  layer,  mixing  of  the  buoyancy  profile 
has  created  an  artificially  strong  gradient  between  the 
level  where  mixing  stopped,  Z„ ,  and  the  next  layer  below,  Z  . 
Thus  salinity  increases  to  compensate  for  the  buoyancy  jump. 
Below  level  Z_,  buoyancy  has  returned  to  a  weak  gradient 
but,  this  level  is  still  in  the  thermocline  of  the  observed 
BT.   Thus  salinity  again  decreases  until  the  temperature 
gradient  becomes  weak.   Below  this,  salinity  again  increases 
because  the  buoyancy  gradient  is  greater  than  ag  times  the 
observed  temperature  gradient.   Thus,  the  second  approach  of 
using  a  typical  salinity  profile  and  adjusting  it  to  the 
bathythermograph  was  judged  superior  in  general.   However, 
notice  this  approach  has  the  disadvantage  of  producing  a 
larger  than  normal  salinity  jump  at  the  base  of  the  mixed 
layer  if  the  typical  salinity  profile  has  to  be  mixed  to  a 
greater  depth  to  correspond  to  the  BT.   Therefore,  at  deeper 
depths,  (>  100  m) ,  it  might  still  be  more  appropriate  to 
use  the  buoyancy  approach,  although  this  suggestion  was  not 
tested  in  this  study. 

The  model  was  initialized  with  a  typical  salinity  profile 
and  bathythermograph  at  OS  PAPA  and  run  for  the  same  13-day 
period  as  above.   The  results  were  cooling  of  the  mixed 
layer  by  0.75   C  and  a  deepening  of  the  mixed  layer  by  5.17  m 
which  agree  relatively  well  with  the  results  obtained  by 
initialization  with  a  BT  alone.   One  might  expect  a  typical 
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profile  to  work  relatively  well  at  this  location  because 
the  overplots  shown  in  Fig.  1  indicate  very  little  variation 
of  salinity. 

By  the  same  token,  one  would  not  expect  a  typical  pro- 
file to  work  as  well  at  OS  HOTEL  given  the  variations  of 
salinity  in  Figs.  2  and  3.   This  was  tested  for  the  same  31- 
day  period  as  above.   Temperature  and  corresponding  derived 
salinity  profiles  are  shown  in  Fig.  2  3.   Notice  that  above 
100  m  the  typical  salinity  profile  has  a  stronger  gradient 
than  the  actual  profile  in  Fig.  15.   Both  the  initial  (Fig. 
23)  and  the  predicted  (Fig.  24)  mixed-layer  salinities  are 
consistently  about  1   /oo  too  fresh.   However,  this  differ- 
ence in  salinity  has  only  a  relatively  small  effect  on  the 
predicted  mixed-layer  temperature  in  Fig.  25,  which  is  simi- 
lar to  the  trend  in  Figs.  17  and  20  with  a  mixed-layer  cool- 
ing of  3.6   C.   Likewise  the  predicted  mixed-layer  depth  in 
Fig.  26  displays  a  similar  deepening  trend  as  in  Fig.  18 
and  is  only  5  m  different  by  the  end  of  the  period  compared 
to  the  model  run  using  the  Nansen  cast  initialization.   Thus 
the  model  has  made  a  good  forecast  even  though  the  typical 
salinity  profile  is  rather  different  then  the  observed  pro- 
file in  the  upper  levels.   The  noticeable  difference  between 
Fig.  26  and  Figs.  18  and  21  is  that  the  mixing,  as  indicated 
by  the  slab  depth  is  strongly  inhibited  by  the  salinity  jump 
at  the  base  of  the  isothermal  depth.   Because  of  the  0.2   C 
temperature  difference  between  the  slab  depths  and  the  mixed- 
layer  depth  definitions,  this  indicates  a  weaker  temperature 
gradient  below  the  well-mixed  layer  than  is  observed. 
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Of  particular  importance  in  Fig.  2 3  is  the  model  pre- 
dicted development  of  a  positive  temperature  gradient  below 
the  mixed  layer.   That  is,  the  model  mixed-layer  tempera- 
ture is  13.78   C  and  the  temperature  increases  to  13.91   C 
at  7.5  m  below  the  isothermal  layer.   A  positive  gradient 
was  observed  in  the  verification  BT.   This  type  of  thermal 
structure  was  also  demonstrated  by  Miller  (1976). 

C.   SOURCES  AND  SINKS  OF  SALINITY 

In  this  section,  the  assumption  of  a  long  term  balance 
between  evaporation  and  precipitation  was  eased  by  including 
surface  salinity  flux  in  the  model.   The  evaporation  rate 
was  derived  from  FNWC ' s  Evaporative  Heat  Flux  (EHF)  field. 
The  precipitation  rate  was  obtained  from  12-hour  accumulated 
precipitation  calculated  by  FNWC ' s  primitive  equation  model 
and  was  assumed  to  be  evenly  distributed  over  the  12  hours. 

Incorporating  the  salinity  flux  at  the  surface,  the 
model  was  initialized  with  the  temperature  and  salinity  pro- 
files in  Fig.  15.   The  resulting  model  profiles,  after  31 
days,  are  shown  in  Fig.  27.   The  significant  changes  in  the 
predicted  salinity  profile  when  compared  to  the  salinity 
profile  predicted  neglecting  (E-P)  are  the  agreement  of  pre- 
dicted and  observed  mixed-layer  salinities  and  the  absence 
of  a  relative  maximum  at  the  base  of  the  mixed  layer  (refer 
to  Figs.  15  and  27).   The  evaporation  rate  minus  the  preci- 
pitation rate  is  shown  in  Fig.  2  8  and  the  corresponding  sur- 
face salinity  in  Fig.  29.   The  surface  salinities  in  Fig.  29 
show  variations  in  response  to  (E-P)  as  opposed  to  the 
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surface  salinities  in  Fig.  16  which  show  only  an  increase 
with  time  due  to  entrainment  of  more  saline  water  at  the 
base  of  the  layer.   Comparison  of  sea-surface  temperatures 
[Fig.  17  without  and  Fig.  30  with  (E-P)]  shows  indistin- 
guishable forecasts.   Thus  sea-surface  temperature  changes 
due  to  inclusion  of  (E-P)  are  apparently  not  significant 
enough  over  a  layer  of  this  depth  for  the  model  to  resolve. 
Finally,  the  slab  depth  responds  differently  under  the  in- 
fluence of  (E-P)  (Fig.  31)  than  when  (E-P)  was  neglected 
(Fig.  18).   The  most  significant  changes  caused  by  (E-P) 
are  the  enhanced  retreats  (decreases  in  slab  depth)  during 
days  3,  10,  and  18  when  precipitation  was  significant,  and 
inhibited  retreats  during  days  8,  13,  17,  22,  25,  and  26 
when  (E-P)  apparently  offset  net  downward  surface  heat  flux. 
Of  particular  note  is  precipitation  on  day  18  which  caused 
a  retreat  that  corresponded  to  the  observations.   The  subse- 
quent slow  deepening  of  the  slab  depth  represents  the  gra- 
dual mixing  of  the  stable  layer  formed  by  precipitation. 

Typical  profile  initialization  runs ,  both  with  and  with- 
out (E-P),  were  made  at  other  ocean  station  locations.   Due 
to  the  sparsity  of  observations,  no  firm  conclusions  can  be 
made  as  to  the  absolute  verification  of  the  forecasts.   How- 
ever, a  relative  comparison  can  be  made  among  results  obtained 
from  the  different  methods  of  initializing  the  model.   These 
results  are  summarized  in  Tables  III  and  IV.  The  length  of 
the  forecast  period  was  determined  by  the  continuous  forcing 
function  data  and  the  longest  period  possible  between  initial 
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and  verifying  thermal  structure  data.   The  different  model 
results  shown  in  Table  III  for  August  indicate  little 
variability  in  predicted  SST  or  MLD  due  to  salinity  effects. 
This  might  be  expected,  in  that  August  is  a  month  of  strong 
downward  heat  flux  and  temperature  effects  would  dominate. 
Thus  one  might  be  cautious  about  including  a  typical  salinity 
profile  under  these  conditions  since  there  is  a  possibility 
of  degrading  the  model  performance.   Included  in  Table  IV 
are  the  results  obtained  during  a  period  of  net  upward  heat 
flux.   One  would  expect  salinity  effects  to  be  more  important 
in  a  deepening  regime.   The  data  in  Table  IV  indicate,  for 
the  most  part,  that  inclusion  of  salinity  results  in  a  warmer 
and  shallower  mixed  layer  than  would  be  the  case  with  den- 
sity as  a  function  of  temperature  alone.   Notice  that,  at 
OS  INDIA  (57N  20W)  and  OS  MIKE  (66N  02E),  the  MLD  exceeded 
300  m  (the  bottom  of  the  model).   This  might  indicate  that 
either  the  "typical"  salinity  profile  used  was  not  really 
representative  of  the  area  or  that  horizontal  advection  can- 
not be  neglected  here.   However,  in  high  latitudes  in  winter 
the  temperature  gradient  is  extremely  weak.   In  this  case, 
the  slab  depth  would  be  more  indicative  of  the  depth  to  which 
the  mixing  processes  were  penetrating.   At  OS  INDIA,  the  pre- 
dicted slab  depths  are  212.5  m  (TYPPRO),  200.0  m  (E-P),  and 
272.0  m  (BT).   At  OS  MIKE,  the  respective  values  are  281.2  m, 
277.5  m,  and  240.0  m.   Although  at  OS  MIKE  the  model  appears 
to  have  run  well  on  temperature  only,  the  initial  BT  was  cor- 
rected for  a  positive  temperature  gradient  which  would  have 
biased  the  potential  energy  of  the  water  column. 
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V.   SUMMARY  AND  CONCLUSIONS 

It  was  the  purpose  of  this  study  to  test  the  effects  of 
salinity  in  a  mixed-layer  model  of  the  ocean  thermal  struc- 
ture.  This  is  desirable  in  that  the  ocean  thermal  structure 
in  many  areas  can  be  explained  only  if  salinity  is  considered 
However,  this  salinity  data  is  not  always  available  even  at 
ocean  weather  stations. 

The  problem  was  approached  using  two  methods:   (1)  ini- 
tialization with  hydrocast  data,  and  (2)  initialization  with 
a  BT  and  a  salinity  profile  derived  from  historical  data 
and  adjusted  to  correspond  to  the  BT.   First,  evaporation 
and  precipitation  were  assumed  to  be  in  balance  over  the 
period  of  integration  and  the  effects  of  salinity  structure 
on  the  thermal  structure  were  observed.   Then,  evaporation 
and  precipitation  were  included  to  observe  the  effects  of 
a  surface  salinity  flux  on  the  thermal  structure.   The  in- 
corporation of  observed  salinity  profiles  into  the  model  did 
not  result  in  significant  changes  in  forecasts  of  SST  and 
MLD  during  the  August  period  tests.   Similarly,  for  the  same 
conditions,  inclusion  of  adjusted  typical  salinity  profiles 
did  not  cause  any  distinct  changes.   However,  the  effects  of 
salinity  structure  were  noticeable  in  deepening  regimes. 
Qualitatively,  model  performance  was  slightly  improved  by 
using  hydrocast  data.   Use  of  adjusted  typical  salinity  pro- 
files and  BTs  also  resulted  in  slight  improvement  even  when 
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the  typical  profile  was  somewhat  different  than  the  observed 
profile.   The  general  effects  of  salinity  structure  were  to 
produce  a  warmer,  shallower  mixed  layer  than  results  when 
density  is  considered  a  function  of  temperature  only. 

Addition  of  the  salinity  flux  due  to  imbalances  of  eva- 
poration and  precipitation  did  not  show  pronounced  effects 
on  SST  or  MLD  in  the  short  term.   However,  a  long  term  up- 
ward buoyancy  flux,  E>P,  tends  to  offset  the  slower  entrain- 
ment  rate  caused  by  the  inclusion  of  salinity  in  the  density 
profile.   The  short  term  effects  of  (E-P)  were  most  evident 
in  the  response  of  the  slab  depth  as  compared  to  the  response 
when  only  the  surface  heat  flux  was  considered.   The  retreat 
of  the  isothermal  layer  is  enhanced  during  precipitation  and 
the  subsequent  deepening  rate  is  slowed.   Additionally,  the 
tendency  of  this  layer  to  shoal  under  downward  heat  flux  may 
be  inhibited  by  evaporation. 

The  addition  of  salinity  effects  to  the  model  is  certain- 
ly desirable.   However,  salinity  data  are  sparse  and  the 
typical  salinity  profiles  must  be  applied  with  an  understand- 
ing of  the  meteorological  and  oceanographical  conditions  of 
the  area.   Thus  one  would  consider  salinity  data  to  be  far 
more  important  in  mid  and  high  latitudes  during  a  deepening 
regime  than  in  low  latitudes  during  strong  heating.   One  must 
also  consider  the  variations  of  salinity  and  the  processes 
causing  these  variations.   In  subpolar  regions,  salinity 
varies  markedly  with  freezing  and  thawing.   In  monsoonal 
areas,  reversal  of  the  wind  and  thus  wind-induced  currents, 
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makes  advective  effects  important.   Further,  an  attempt  to 
apply  the  model  ocean-wide  would  require  salinity  data  in 
many  areas  where  it  is  not  now  available. 
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APPENDIX  A 

HYDROCLIMATOLOGICAL  DATA 
RETRIEVAL  PROGRAM 


The  Hydroclimatological  Data  Retrieval  (HYDAT)  program 
was  designed  by  Ocean  Data  Systems,  Inc.  (1977)  to  access, 
process,  display,  and  store  hydroclimatological  data  bases. 
Among  the  outputs  available  from  HYDAT  is  an  option  that 
produces  a  typical  temperature  profile.   This  typical  profile 
is  obtained  by  specifying  the  month  (or  months)  and  the  loca- 
tion of  interest.   The  location  specified  is  taken  to  be  the 
midpoint  of  a  circular  or  square  search  area,  with  the 
original  search  distance  from  the  midpoint  determined  by  the 
user.   If  at  least  six  reports  are  not  found  in  this  area, 
the  area  is  expanded  by  fixed  increments  until  six  or  more 
reports  are  found  or  until  the  maximum  search  distance, 
either  default  or  specified,  is  reached. 

If  the  search  is  successful,  each  report  found  is  sepa- 
rated into  three  parameters:  (1)  the  sea-surface  temperature, 
(2)  the  deep  temperature,  a  temperature  at  an  arbitrary  depth 
(~2  00  m) ,  and  (3)  the  heat  content  as  computed  from 

N 


>   [T-(^)]'AZ  (A-l) 


i  =  l 

where  D  is  the  deepest  level,  N  is  the  number  of  layers  to 
depth,  D,AZ  is  the  thickness  of  each  layer,  and  T  is  the 
mean  temperature  in  the  layer.   Any  sounding  not  extending 
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to  D  is  not  considered  in  the  selection  of  the  typical  pro- 
file.  The  values  of  each -parameter  are  then  screened  to 
remove  isolated  extreme  values  and  the  frequency  distribu- 
tion is  checked  for  the  number  of  modes  (<_  3)  present.   All 
three  parameters  must  have  the  same  number  of  modes.   HYDAT 
then  computes  and  prints  the  median  and  first  and  third 
quartile  values  for  each  mode  of  each  parameter.   The 
typical  profile  in  each  mode  is  determined  by  comparing  the 
parameter  values  of  each  profile  with  the  median  values  and 
determining  a  deviation,  D., 

W. | P.  -  M. | 
D.  =  -! ijj- L_  ,  1=1,2,3       (A-2) 


where  W.  is  a  weighting  factor  assigned  to  the  individual 
parameter,  P.  is  the  value  of  the  parameter,  and  M.  is  the 
median  value  of  the  parameter  within  any  one  mode.   The 
typical  profile  is  then  the  report  with  the  lowest  sum  of 
deviations : 


Sum  =  Z  D.  (A- 3) 
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TABLE    II.       Comparison    of    STP    casts  with    BTs    at    OS    PAPA, 

November-December    1974.  Observations    are 

within    one   hour   and   ten  nautical   miles    of 
each   other. 


b^m 


ssrfcrp) 


ZSMcV 

8.11  °c 

£4ICW 

7.40 

'Urn    0<W 

7,87 

2?A>o/ 

7.^1 

01  bee. 

7.84 

6l&:c 

7.13 

o3D«?c 

7.(d1 

C^Qeu 

7.8| 

OtOeo 

74,8 

55T16D 


fAoxn  <^.V\e^<iv\ce.<> 


8.D°e 
7.<o 

7.6, 
7.(o 

7  7 

7.4 

7,5 

1M 


LT 


6.11  C 
0.3 

0.2-7 
O.L1 


/-u-DlsrP) 


6o.  J     m 
7©.^* 
7Z..I 
11.0 

(cS,  I 

If.  3 
7b.  (d 


O.ZC°£ 


MM)(feT) 


7k.  2.   »n 
73.8 
6b- 7 

61.3 
71.0 

$6.3* 


2iD 


-10.7  m 
-Z.cl 

-12.4 

o.3 

-jo.'Z. 


-i/.4*\ 


Xnd-ioxieji  pc6j"/iv/e  cif&cite>o4  «^tcu-vi4 ef ed  c\i6ove  MU), 


41 


TABLE  III. 


Comparison  of  forecast  SST  and  MLD  among 
different  model  initialization  runs  and 
with  observations  for  August  1974.  HC  is 
hydrocast  data,  TYPPRO  is  BT  and  typical 
salinity  data,  (E-P)  indicates  inclusion 
of  surface  salinity  flux,  and  BT  is  tem- 
perature data.   Verifying  data  are  from  BTs 
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TABLE  IV.   Comparison  of  forecast  SST  and  MLD  among 
different  model  initialization  runs  and 
with  observations  for  November-December  19  74 
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Figure    1.       Overplot    of  hydrocast    data    for    OS    PAPA    for 
the    period    21    November-31    December    1974. 
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Figure    2.       Overplot    of   hydrocast    data    for    OS    HOTEL    for 
the   period    21   November-31    December    1974. 
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Figure    3.       Overplot    of   hydrocast    data    for   OS    HOTEL    for 
the    period    8-28    August    1974. 
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Figure  10.   Depth  vs  temperature  at  OS  PAPA  -  18  GMT  2  3 

November  and  18  GMT  06  December  1974,  BT  data. 
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Figure    15. 


Depth    vs    temperature    and    salinity    at    OS    HOTEL 
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1974,    hydrocast    data,    neglecting    (E-P). 


58 


IS) 


-«» • 


_, 4-j 


■«•■- ««f>  - 


!       J      | 

O      CQ 

•i -s- p» 


I       1  i 


4 + — t  -- 


..+-..!.. 


csa 


(S) 
CM 


C2 

CO 
C-J 


IS 

■ 

C--J 


"S3 

■ 

C--J 
CM 


CS) 


CS) 
C-.J 


IS) 
CO 

IS) 

to 

1 — 1 

IS) 

— I 

S3 

i 

cm 


eg 


... 


;v I 


C3 
03 

CS2 


IS) 
CS) 

cm 

"S3 

cs 


csa 

CSl 

CD 

ro 


IS) 

en 

CO 


Q 
CO 

on 


IS) 
LD 

oo 


C3 
CO 

m 
en 


en 


is 

li~J 

ro 


IS) 


on 


csa 

CXI 
Lf) 

no 


cs 


LD 

m 


S3 

cs 
ro 


-*>«■ 
^ 


O 
M 

CO 
CM 


O 
CC 
L~ 

CO 

>- 

cr 

Q 


o 

-p 

,Q 
6 

> 

O 

s    • 

OO    Qh 

CM       I 

W 

E-»  ~ 
S 
CD    bO 

C 
J-  -H 
rH    +-> 

O 
I      CU 


W  d) 

H  C 
O 

EC  « 

co  r- 
O  en 

rH 

rd    P, 
CI) 

■H  £ 
•H    d) 

c  o 

•H  CU 
-H  Q 

rd 

co  en 

CM 

0) 

O  Eh 
<3  S 
4h  CD 
P. 

CO  rH 


<X> 


CU 

bO 
•H 


Cldd)  A1INI1US 


59 


o 


M 

CO 
C\J 


o 
cc 
Li- 
en 

>— 
cr 

Q 


CO 

1 

CM 

O 
P. 

H  T3 

S 

>i 

CD     rC 

J" 

s 

1 

o 

14H 

-J 

J- 

W 

r-~ 

Eh 

CX> 

O 

rH 

E 

P, 

00 

CD 

o 

6 

-H 

dJ 

rd 

o 

CD 

CD 

O 

P. 

3 

en 

•H 

CM 

• 

rcj 

c 

fc  ^ 

0 

0)  S 

•H 

£XO 

■H 

e 

ft! 

a) 

J- 

N 

+-> 

rH 

•H 

in 

0 

rti 

a; 

-P 

•H 

>> 

-H 

d 

P. 

•H 

!-H 

CD 

c 

1 

^ 

■H 

TD 

6 

CD 

cd 

+-> 

X 

> 

CO 

•H 

0 

nj 

s 

^ 

a 

CO 
Li 

bO 

•H 
U-, 


(3  030)   dW31 


60 


eg 
en 


CO 


CO 

-•J 


CM 


N 

Xj 

C^ 

C3 

^J, 

IS 

o 

C\l 

.<£. 

C3 

hsJ 

• 

-*t 

QO 

* — i 

CO 

CM 

13 

CO 

ZI 

— ' 

o 

cc 

q 

Ll. 

^ 

co 

_>— 

CE 

"S3 

Q 

CM 

— 1 

IS) 

S2 

0 

c 

■p 

o 

•H 

h 

-H 

0) 

ft) 

,q 

N 

i 

•H 

QJ 

H 

> 

rd 

O 

•H 

s 

•H 

•H 

CO 

C 

CM 

•H 

E-i 

+J 

£ 

W 

CD 

rd 

a 

J- 

o 

rH 

U 

TD 

1 

>, 

rC 

-J 

w 

6 

H 

0 

o 

& 

EC  4h 

00 

Zi- 

o 

eh- 

en 

+-1 

H 

id 

P. 

A 

(D 

+-> 

fl 

P.  6 

0) 

Q) 

T3 

O 

CD 

& 

O 

Q) 

>^ 

en 

rd 

CM 

H 

1 

H 

TD  2 

0) 

CD 

X 

•H 

J- 

S  rH 

CD 
Sh 

bO 
•H 
Pn 


IS) 

CSJ 

cs 

(Si 

>SI 

cs 

EJ 

• 

• 

• 

• 

• 

. 

• 

LD 

CO 

c*» 

CO 

en 

E3 

-s* 

T— t 

C\J 

ro 

^ 

m 

f^ 

(W)   Hld3Q 


61 


20. 0 


CD 

UJ 
Q 


INITIAL 


15.0 


-I  -    -    MODEL 
I-    —    OBS  | 


10.0 


5.0 


0.8 


100.0 


150.0 


300. 0 


DEPTH   CM) 


Figure    19.       Depth    vs    temperature    at    OS    HOTEL   -    14    GMT    2  8 
November   and   14    GMT    29    December,    BT   data. 
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Figure  22.   Schematic  of  a  derived  salinity  profile 
assuming  a  typical  buoyancy. 
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Figure    2  3. 


Depth    vs    temperature    and   salinity    at    OS    HOTEL 
-    14    GMT   2  8    November   and    14    GMT   2  9    December 
1974,    BT   and   typical    salinity    data,    neglecting 
CE-P). 
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Figure    2  7.       Depth   vs    temperature ' and    salinity    at    OS    HOTEL 
-    14    GMT   2  8    November   and    14    GMT    2  9    December 
19  74,    hydrocast    data. 


70 


o 


M 

■< — i 

oo 

CM 


O 
CC 
U_ 

CO 

en 

Q 


CD 


O 
-H 

Sh 

QJ 
^1 
S 

0) 

> 

O 

oo 

CM 

Eh 

CD 


I 

•J  • 

W  J- 

H  r» 

O  <T> 

K  rH 

CO  in 

O  (U 

,Q 

+J  6 

tti  QJ 
0 

^  QJ 

0<  Q 

W  cn 

v_/  CM 


<X) 
CM 

0) 

& 

bO 
•H 


0T«yH/WG(d-3) 


71 


4—f 


•+-U 

W  ! 
Q  I  00 
■•OT'CQ- 
S  i  O 


i  + 


..+..;.. 


13 

csa 


CO 


13 

■ 

CO 

r-vj 


13 


is 

O.J 

r\J 


CSJ 


!3 
CD 

3 
CO 


CO 


CN 

csa 


ea 

CO 

no 


is 

cn 
IT) 


F3 
CO 

LO 
CO 


IS 

cn 
oo 


(3 


cn 

(TO 


IS 

cn 
cn 

CO 


4"0 


G3 

cn 

cn 
00 


ca 

CSI 

cn 
ro 


S 


cn 
ro 


O 

CSl 

cn 
no 


> 

O 

rvi 
-** 

•< — i 

00 

CM 


o 

■H 

Sh 
cu 

■i 

cu 
> 

o 

cn 
cnj 

Eh 
CD 


—4 

ro 

J 

or 

w 

13 

U_ 

Eh 
O 

PC     • 

~W 

CO 

zr 

>- 

cn  r- 

3 

cx 

Q 

O  cn 

rH 
-M 

cn 

<tf    Sh 

x— i 

cs 

•H    0) 

1 

c  o 

«— 1 

sali 
9    De 

13 

CM 

CU 

CO 

O   H 

M-t  CD 

CO    rH 


cn 

OJ 


cu 

Sh 

bO 
•H 
Cm 


(ldd)   AlINIlbS 


72 


f—\ 


CO 
CM 


Q 

or 
u_ 

en 
>- 

cr 
a 


u 

0) 

> 

o 

CO 

Eh 

CD 


W 

o 
o 

■H 

fti 

h 

■H 
it) 
U 

CD 

a  o 
6  a) 
<u  p 

+-» 

en 
U   CM 
<L> 
>,Eh 

id  s 

H  CD 


i> 

en 


CO 

£ 

CD 


-a 

■H 


O 


o 

CO 

CD 
Sh 

=3 

bO 

•H 

p4 


(3  930)   dW3i 


73 


I    iJ         i 
»..«!  — . j «  _;  „ j,  \ _., .pa f  ~w.~ j.._- 

ji/  !        <      I    Q  tO 

t1{             i       -J     i    O         PQ 
... ...  Tj  ,..._..u:j.„._....s;.. n....... 

i 
I,' 


C3 
00 


sa 


r-0 


C3 

0 

i 

+j 

■«% 

<M 

0 

(N 

-■vf 

0) 

> 

l>- 

0 

rs 

^^ 

^> 

CO 

ca 

r    1 

r^ 

2H 

Csl 

Eh 

en 

M 

S 

xj 

CD 

CO 

! 

zf 

00 

H 

C\J 

(9 

1 

rjO 

5Z 

w 

■^-4 

O 

tr 

Fh 

■3 

Li- 

O 

-*t 

en 

• 

T— 1 

%_ 

CO  j- 

/r 

O  c> 

CE 

CD 

IS 

Q 

■M   rH 

,— 4 

rd 

cs 

ft  e 

• 

(D    CD 

is 

T3    a 

a) 

&  Q 

a) 

O 

ft!  cn 

00 

ed-1 

GMT 

C3 
■ 

•H  zr 

f£) 

S  H 

00 

CD 

bO 
•H 
Ph 


CS1 

C3 

s 

• 

• 

• 

r-. 

CO 

CO 

«— 1 

<M 

(WJ    Hld3Q 


74 


LIST  OF  REFERENCES 


Camp,  N.  T. ,  1976:   Role  of  strong  atmospheric  forcing 

events  in  the  modification  of  the  upper  thermal  struc- 
ture during  the  cooling  season.   Ph.D.  thesis,  Naval 
Postgraduate  School,  175  pp. 

Camp,  N.T.  and  R.  L.  Elsberry ,  1978:   Ocean  thermal  response 
to  strong  atmospheric' forcing.  II.  The  role  of  one  di- 
mensional processes.   J.  Phys.  Oceanogr. ,  To  appear  in 
a  forthcoming  issue. 

Elsberry,  R.  L. ,  T.  S.  Fraim,  R.  N.  Trapnell,  1976:   A 
mixed-layer  model  of  the  oceanic  thermal  response  to 
hurricanes.   J.  Geophys .  Res . ,  81 ,  1153-1162. 

Elsberry,  R.  L.  and  N.T.  Camp,  1978:   Ocean  thermal  response 
to  strong  atmospheric  forcing.  I.  Characteristics  of 
forcing  events.   J.  Phys.  Oceanogr.,  To  appear  in  a 
forthcoming  issue. 

Garwood,  R.  W. ,  Jr.,  1977:   An  oceanic  mixed-layer  model 

capable  of  simulating  cyclical  states.   J.  Phys.  Oceanogr. 
7_,    455-468. 

Haney ,  R.  L.  and  R.  W.  Davis,  1976:   The  role  of  surface 
mixing  in  the  seasonal  variation  of  the  ocean  thermal 
structure.   J.  Phys.  Oceanogr.,  6_,  504-510. 

Institute  of  Ocean  Sciences,  1975:   Oceanographic  Observa- 
tions at  Ocean  Station  P.   Unpublished  manuscript. 

Johnson,  W.  F. ,  1977:   Upper  ocean  thermal  structure  forecast 
evaluation  of  a  model  using  synoptic  data.   M.S.  thesis, 
Naval  Postgraduate  School,  47  pp. 

Kraus ,  E.  B.  and  J.  S.  Turner,  1967:   A  one-dimensional  model 
of  the  seasonal  thermocline.  II.  The  general  theory  and 
its  consequences.   Tellus ,  19 ,  98-106. 

Miller,  J.  R. ,  1976:   The  salinity  effect  in  a  mixed  layer 
ocean  model.   J.  Phys.  Oceanogr.,  6,  29-35. 

Ocean  Data  Systems,  Inc.,  1977:   Hydroclimatological  Data 
Retrieval  Program.   Unpublished  manuscript. 

Thompson,  R.O.R.Y.,  1976:   Climatological  numerical  models 
of  the  surface  mixed  layers  of  the  ocean.   J.  Phys. 
Oceanogr. ,  6,  496-503. 


75 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California   9  3940 

3.  Director  Naval  Oceanography  and  Meteorology      1 
National  Space  Technology  Laboratories 

Bay  St.  Louis,  Mississippi   39520 

4.  Professor  R.  L.  Elsberry,  Code  63Es  8 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California   9  3940 

5.  Professor  G.  J.  Haltiner,  Code  63Ha  1 
Chairman,  Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California   93940 

6.  Department  of  Oceanography,  Code  6  8  1 
Naval  Postgraduate-  School 

Monterey,  California   93940 

7.  Lieutenant  R.  A.  Paulus  2 
Naval  Ocean  System  Center,  Code  035 

271  Catalina  Blvd. 

San  Diego,  California   92152 

8.  Commanding  Officer  2 
Fleet  Numerical  Weather  Central 

Monterey,  California   9  3940 

9.  Naval  Oceanographic  Research  and  1 
Development  Agency,  Code  0201 

Bay  St.  Louis,  Mississippi   39520 

10.  Professor  R.  L.  Haney ,  Code  6  3Hy  1 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California   9  3940 

11.  Professor  R.  W.  Garwood,  Code  68Gd  1 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California   93940 

76 


12.  Lieutenant  V.  F.  Looft 
Naval  Weather  Service 
Environmental  Detachment 
Moffett  Field,  California   94035 

13.  Lieutenant  S.  M.  Schrobo 
OA  Division 

USS  KITTY  HAWK  (CV-6  3) 
FPO  San  Francisco   96601 

14.  Lieutenant  M.  G.  Salinas 
Naval  Support  Forces  Antarctica 
Met.  Div. ,  Code  030 

FPO  San  Francisco   96601 


77 


175439 


Tresis 
P2-744 


Thesis 

P2744 

c.l 


H33 


■  -  an 

■*  7 

Paul  us  " 

Salfnfty  effects    fn 
an   ocean  c  miw»w    1 
mode?.  m'*e<Mayer 

31  jUL  2  6  4  9  3 

3  3  9  0  B 


Thesis 

P274H 
c.l 


Paul  us 

Salinity  effects  in 
an  oceanic  mixed-layer 
mode  I . 


?.  75W9 


thesP2744 


SmZmm,S  m  an  0Cean,c  mixed-lav 


3  2768  001  98094  9 

DUDLEY  KNOX  LIBRARY 


